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Abstract 

We show that the N — 2 superextended ID quantum Dirac delta potential problem is char- 
acterized by the hidden nonlinear su(2|2) superunitary symmetry. The unexpected feature of 
this simple supersymmetric system is that it admits three different Z2-gradings, which produce 
a separation of 16 integrals of motion into three different sets of 8 bosonic and 8 fermionic oper- 
ators. These three different graded sets of integrals generate two different nonlinear, deformed 
forms of su(2|2), in which the Hamiltonian plays a role of a multiplicative central charge. On the 
ground state, the nonlinear superalgebra is reduced to the two distinct 2D Euclidean analogs of 
Oh' a superextended Poincare algebra used earlier in the literature for investigation of spontaneous 

supersymmetry breaking. We indicate that the observed exotic supersymmetric structure with 
three different Z2-gradings can be useful for the search of hidden symmetries in some other 
quantum systems, in particular, related to the Lame equation. 
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The Dirac delta potential plays a prominent role in modeling diverse physical systems and phenom- 
■ ena [1, 2, 3, 4]. Recently, it was observed that the simplest bosonic ID quantum Dirac potential 

problem is characterized by a hidden nonlocal N = 2 supersymmetry in which a reflection plays 
a role of the grading operator [5]. The corresponding supersymmetry is exact in the case of the 
attractive potential and is spontaneously broken in the repulsive case. 
^ . A priori it is clear that if the system is superextended by introduction of the discrete (spin) 

degrees of freedom, its explicit supersymmetry should enlarge. The purpose of the present paper 
is to investigate how the explicit N = 2 supersymmetry of the system unifying the attractive and 
repulsive potential cases is enlarged by the nonlocal hidden supersymmetry. 

We show that the N = 2 superextended ID quantum Dirac potential problem is characterized by 
a nonlinear su(2|2) superunitary symmetry. The unexpected feature of this simple supersymmetric 
system is that it admits three different possible choices of the grading operator T, a usual one given 
by T = o"3, and two additional given by T = R and T = R0-3, where R is a reflection operator. 
These three choices provide us with three different separations of the total set of 16 integrals of 
motion into 8 bosonic and 8 fermionic operators. The nonlinearity of the superalgebra formed by 
the integrals of motion is similar to the nonlinearity of the symmetry algebra appearing in the 
Kepler problem and associated there with a hidden symmetry provided by the Laplace-Runge-Lenz 
vector [6]. We find that the three choices for the grading operator give rise to the two different 
forms of nonlinear, deformed su(2|2) supersymmetry, with the Hamiltonian playing the role of the 
multiplicative central charge of the superalgebra. It is worth to note here that for T = R0-3, the 
usual supercharges (2.3) play a role of the bosonic generators of the superalgebra, whose form in 
this case is the same as for the choice T = 0-3 which identifies (2.3) as fermionic integrals. 
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On the ground state, the two types of superalgebras are reduced to the 2D Euclidean versions 
of the two essentially different cases of a nonstandard superextension of the Poincare algebra, which 
was used earlier by Gershun and Tkach in investigation of spontaneous supersymmetry breaking 
[7]- 

The paper is organized as follows. In Section 2 we identify the set of integrals of motion of the 
system, and the form of the nonlinear superalgebra generated by them depending on the choice of 
the grading operator. In Section 3 we consider the action of the integrals of motion on the energy 
eigenstates. In Section 4 we identify the nature of the nonlinear supersymmetry of the system. 
In Section 5 we investigate its reduction on the ground state. Finally, section 6 is devoted to the 
discussion and concluding remarks. 



2 Integrals of motion and grading operators 

Consider an N = 2 supersymmetric extension of the ID quantum Dirac potential problem described 
by the Hamiltonian 

H = p 2 + 1 + 2(35{x)a 3 . (2.1) 

We choose the units h = c = 1, and put the mass of the particle m = ^. In this units, coordinate, 
momentum, energy and the real parameter (3 > are dimensionless. The system is constituted 
by the two supersymmetric partners corresponding to the attractive (lower component) and the 
repulsive (upper component) potentials. A usual, the explicit supersymmetry of the system (2.1), 

{Q a ,Q b } = 2H5 ab , [Q a ,H}=0, a = 1,2, (2.2) 

is generated by the supercharges 

Qi = P&1 + f3e(x)a 2 , Q2 = i<??,Qi- (2.3) 

Here e(x) is the sign function, e(x) = +1 (—1) for x > (x < 0) and e(0) = 0, -^s.{x) = 25(x), 
which can be treated, e.g., as a limit e(x) = lim^oo tanh Xx, and then 5(x) = lim^oo A/ cosh 2 Xx. 
The identities for the Pauli matrices, GjOk = Sjk + iejkici, give an equivalent representation for 
(2.3), 

Qi = <?i {p + i(3e(x)a 3 ) , Q 2 = -cr 2 (p + i(3e(x)a 3 ) . (2.4) 
For the N = 2 superalgebra (2.2), a ^-grading operator T, 

r 2 = l, [T,H}=0, {T,Q a } = 0, (2.5) 

is identified usually with the diagonal a-matrix, 

r = (73. (2.6) 

This choice is not a unique, however, and the identification of the grading operator with the 
reflection, 

T = R, (2.7) 

Ripix) = ip(—x), is also consistent here with relations (2.5). 

In accordance with [5], the system has also the integrals of motion 

Qi=p + i(3e(x)Ra 3 , Q 2 = iRQi, (2.8) 

which, like R, are nonlocal operators. Taking into account the identity 5{x)R = S(x), one finds that 
like usual supercharges (2.3), they are square root of the Hamiltonian, Q\ = Q\ = H. Integrals 
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(2.8) are the Z2-odd operators (supercharges) for the choice (2.7), but have to be treated as Z2-even 
operators for (2.6). In the first case, they generate another copy of the N = 2 superalgebra. 
There is a third possibility for identification of the grading operator, 

T = Ra 3 . (2.9) 

With respect to this grading, operators Q a are odd, while Q a are even. 

Let us see now what is the whole set of even and odd integrals of motion of the system for each 
of the three different grading operators we have just specified. 

2.1 Grading T = R 

Since with respect to (2.7) both sets of integrals (2.3) and (2.8) are odd operators, let us fix this 
grading and analyze further the symmetries of the system. With respect to the other two gradings, 
one of the two sets of integrals (2.3) and (2.8) should be treated as odd, and the other as even. 
These alternative choices of the grading will be considered separately, but we shall see finally that 
in all the three cases the structure of the resulting complete supersymmetry is the same, modulo 
the Hamiltonian. 

For the grading (2.7), the anticommutators between the integrals (2.3) and (2.8) should be 
calculated. This gives 

{Q a ,Qi} = 2S a , {Q„,Q 2 } = 0, a = 1,2, (2.10) 

where 

S 1 =<r 1 H-Pe(x)Q 2 (l + R), S 2 = ia 3 S 1 . (2.11) 

Hermitian operators (2.11) should be treated as new even integrals of motion. Then, the process 
has to be continued: it is necessary to calculate the commutators of these even integrals of motion 
between themselves and with the odd integrals, etc. As a result, we obtain the complete list of odd 
(Fi, . . . , Fg) and even (H, R, Si, £2, B±, . . . B4) integrals of motion, which can be represented in 
terms of the basic 'building block operators' Q\, Qi, Si, a 3 and R. They are shown in Table 1. 



Fermionic 
integrals 


Fi = Qi 


F 2 = Ra 3 Qi 


F 3 = Q 2 = icr 3 Q 1 


F 4 = iRQi 


F 5 = iRQi = Q 2 


F 6 = iRa 3 Q 1 


F7 = 


F 8 = Qi 


Bosonic 
integrals 


H 


t = r 


Si = 03 


£2 = R(T 3 


Bi = Q1Q1 = Si 


B 2 = icr 3 Si = S 2 


B 3 = RSi 


B4 = iRa 3 Si 



Table 1: Integrals of motion, grading V = R 

The anticommutation relations between the fermionic integrals are presented in Table 2. 

The commutation relations between bosonic integrals and between bosonic and fermionic oper- 
ators are listed in Tables 3 and 4. 

Here we have introduced the parameter A to include also the other two cases of ^-grading. 
In the present case of T = R, A = 0. Due to the special choice of the indices in the definition of 
fermionic operators F and bosonic operators B, the tables of (anti) commutation relations have a 
simple block form. The Hamiltonian H commutes with all the fermionic and bosonic operators. 
The grading operator T = R commutes with all the bosonic integrals. 

2.2 Grading T = cr 3 

For the choice of the ^-grading operator T = a 3 , as we noted, the integrals (2.3) are, again, the 
fermionic operators, but the integrals (2.8) are even operators. As a result, now the operators S a 
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Table 2: Fermion-fermion anticommutation relations 





Fi 


F 2 


^3 


F 4 


^5 


F e 


F 7 


F 8 


r 


-2iF A 


-2iF 3 


2iF 2 


2iF x 


2iF 8 


2iF 7 


-2iF 6 


-2iF 5 


Si 


-2iF 3 


-2iF 4 


2iF x 


2iF 2 
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2iF 7 


2iF 8 


-2iF 5 


-2iF 6 


Bi 





2iF 6 H 1 - x 


-2iF 7 H 1 ~ x 








-2iF 2 H 


2iF 3 H 





B 2 


2iF 7 H 1 - x 








2iF 6 H 1 - x 





-2iF 4 H 


-2iF 1 H 





B 3 


-2iF b H l ~ x 








2iF 8 H 1 - x 


2iF 1 H 








-2iF 4 H 


B 4 





-2iF 8 iJ 1 - A 


-2iF b H 1 ~ x 





2iF 3 H 








2iF 2 H 



Table 3: Boson- fermion commutation relations 



are odd operators appearing from the commutation relations of (2.3) with (2.8). The identification 
of the fermionic and bosonic integrals of motion for this case is presented in Table 5. 

These operators satisfy the (anti) commutation relations of the same form as in the previous 
case, but now with A = 1. 

2.3 Grading T = Ra 3 

The choice T = Ra 3 is similar to the choice (2.6). The difference with the previous case is that now 
the integrals (2.3) are identified as even operators, while the integrals (2.8) have a nature of odd 
operators. The integrals S a , together with two other integrals related to them, play again here the 
role of odd operators. The identification of the fermionic and bosonic operators is given in Table 6. 

These operators satisfy the (anti) commutation relations exactly of the same form as for grading 
(2.6), i.e. the parameter A is assigned here the value A = 1. 

2.4 Identification of operators 

The identification of the fermionic and bosonic operators has been realized by us in a special way 
which guarantees the same form (modulo H) of the (anti)commutation relations between them for 
the three possible choices of the Z^-grading operator. The identification we use corresponds to the 
following simple procedure. 

Step 1. Choose the grading operator T from the set of the three Hermitian operators 

{R, a 3 , Ra 3 }. (2.12) 
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2iY ll H 2 ~ x 






Table 4: Boson-boson commutation relations 



Fermionic 
integrals 


Fi = Qi 


F 2 = -RasQi 


F 3 = -iRQi 


F 4 = Q 2 = ia 3 Q\ 


F 5 = RS! 


F 6 = -S! 


F 7 = ia 3 Si 


F 8 = -iRa 3 Si 


Bosonic 
integrals 


H 


Si = -R 


T = a 3 


S2 = —Rcr 3 


Bi = -iRa s Qi 


B 2 = -cr 3 Qi 


B 3 = -iRQx 


B 4 = —Qi 



Table 5: Integrals of motion, grading T = a 3 



Fermionic 
integrals 


Fi = Qi 


F 2 = —a 3 Qi 


F 3 = -iRQi 


F 4 = iRo 3 Q\ 


F5 = RSi 


F 6 = -5i 


F-j = iRa 3 S\ 


F s = -ia 3 Si 


Bosonic 
integrals 


H 


Si = -R 


S2 = — CT 3 


T = Ra 3 


B\ = -ia 3 Qi 


B 2 = -Ra 3 Qi 


B 3 = -iRQi 


B 4 = —Qi 



Table 6: Integrals of motion, grading T = Ra 3 

Step 2. Select any Hermitian fermionic operator F±, {Fi,T} = 0, from the set of the eight odd 
integrals with the following properties: {Fi,Fi} = 2H and [£2,-^1] = 0. Here we denote by 
S2 / T a bosonic operator from the set (2.12) which commutes with F\. The third operator 
from (2.12) is denoted by Si. The operators Si and S2 are defined up to a sign. 

Step 3. With the commutation relations from Table 3 of F\ with T and Si, the fermionic operators 
F 2 , F 3 and F 4 are fixed. 

Step 4. Choose any other fermionic operator F5 that satisfies {F±, F5} = and {F5, -F5} = 2H l+x . 

Step 5. Repeat Step 3, but changing F\ — ► F5 and Si — ► S2 to obtain Fq, F-j and F§. 

Step 6. With the anticonmutators of the Tables 2 we identify the bosonic operators Bi, B 2 , B 3 
and B 4 . 

3 The action of the integrals of motion on the energy eigenstates 

We have three types of the 'basic' integrals of motion, Q a , Q a , S a , in terms of which other integrals 
are obtained by multiplication with the even integrals R and a 3 . Here we will trace out the difference 
between the nature of these operators by applying them to the physical states which we choose 
finally to be the eigenfunctions of the operators H, R and a 3 . 
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In the repulsive case (the upper component subsystem), the Hamiltonian is H = p 2 +2/35(x)+f3 2 , 
and the eigenfunctions with energy E = k 2 + /3 2 > /3 2 , k > 0, are given by 

4,+ (x) = (e ikx + re- tkx )e(-x)+te tkx Q(x), ^-(x)=^(-x), (3.1) 

where 0(x) is the Heaviside step function. These wavefunctions correspond to scattering states 
with the plane waves incoming, respectively, from — oo and +00, with reflection and transmission 
coefficients given by [5] 

Let us construct from these states the eigenfunctions of a fixed parity, 

ip$.(x) = ^ ~ ® Wt + ^k) = k cos kx + P< x ) sin kx > ( 3 - 3 ) 

^- {x) = " ^ = s[nkx > (3 - 4) 

The attractive (lower component) case is obtained by the change (3 — > — /? in all the previous 
formulae, H = p 2 - 2/35(x) + /3 2 , 

+0*0 = & cos ^ X — P £ ( x ) ^ 1{ X ) = sinkx. (3-5) 

In addition to the scattering states, this subsystem has a unique bound state described by a nor- 
malized wave function, 

iP { -\x) = y/pe-PW. (3.6) 

Due to a special value of the constant term in the Hamiltonian, the energy of this bound state is 
equal to zero. 

Summarizing, for the system (2.1) we have the scattering eigenstates, 

*£>(*) = ( ^te + PeW^hx ^ ^ {x) = ^ sinkx ^ (3 7) 

, Ir,- h V ^(x) = V^TW ( °,), (3.8) 
V k cos kx — p£\x) sin kx J \ suite y 



and the bound state of zero energy, 



*» W =U'-*' )' <3 ' 9) 

We have introduced here a convenient normalization to simplify the relations that follow. 

The ordinary supercharges (2.3) change the parity and the upper-lower subspaces of the eigen- 
functions, 

= ±i^k 2 + p 2 ^\ Qi^ ] = TWk 2 + P 2 *±\ (3-10) 

Q 2 ^ = ^WTW^\ Q 2 ^ ] = -^JWTW^±\ (3-11) 

where for simplicity we have omitted the index k. The supercharges of the hidden supersymmetry 
(2.8) change only the parity, 

= ±i^/WTW^\ Q 2 *L Q) = y/tf + prf" \ « = +,-, (3.12) 
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The S a operators change only the upper-lower subspaces of eigenfunctions, 

SitfW = (k 2 + (5 2 )^\ S 2 *W = ±i(fc 2 + /3 2 )^ T) , « = +,-. (3.13) 

The unique bound state is annihilated by the integrals of motion Q a , Q a and S a , 

Q a *S _) =0, Qa^ } = 0, S a tf< _) =0. (3.14) 

As a consequence, it is annihilated by all the odd integrals of motion Fi,...,F& and all the even 
integrals B\,...,B± (as well as by H). 

4 Identification of supersymmetry, generic case 

To clarify the nature of the supersymmetry generated by the complete set of eight fermionic and 
eight bosonic integrals of motion, we define the following linear combinations of the even operators, 

# } = ^(Si± ^ 2 ) = l^iU ± , (4.2) 
where n± = |(1 ± V) are the projectors. These operators satisfy the relations 

= i4 ±] H 2 -\ (4.3) 
= ie ab P^\ a, b = 1,2, (4.4) 

= lA^iP^] = [jt\H +) ] = [Jt\jt ] ] = 0. (4.5) 

The commutation relations (4.3), (4.4) correspond to a direct sum of two deformed su(2) al- 
gebras. In particular, relations (4.3) are reminiscent of the commutation relations satisfied by 
the components of the Laplace-Rungc-Lenz vector in the quantum Kepler problem. Being re- 
duced to the energy eigensubspace of nonzero eigenvalue E > (3 2 > 0, the rescaled operators 
ji ± - ) = P^ /E 1- ^ and generate the Lie algebra su(2) © su(2) and satisfy the relations 

= |n + , j| ^ j\ ^ = |n_, where the summation in i = 1, 2, 3 is assumed. These relations 
mean that the two common eigenstates of the Hamiltonian with energy E > (3 2 > and of the 
grading operator T with eigenvalue +1 carry spin-1/2 representation for J- + \ and are the spin-0 
states for . The states with T = — 1 carry spin one-half representation for j\ ~\ and are of 

spin zero for . Fermionic integrals mutually transform the states from these eigenspaces of the 
grading operator. In accordance with the total number of independent fermionic generators, the 
energy subspace with E > (3 2 > carries an irreducible representation of the sit(2|2) superunitary 
symmetry, which is a super symmetric extension of the bosonic symmetry u(l) © su(2) © su(2), 
where the u(l) subalgebra is generated by the grading operator, see Ref. [8]. Having in mind that 
the Hamiltonian appears in a generic form of the superalgebra as a multiplicative central charge, 
we conclude that the system possesses the nonlinear su(2|2) superunitary symmetry in the sense 
of Refs. [6, 9, 10]. 

In order to understand better the identified symmetry at the algebraic level, let us define the 
following linear complex combinations of the even operators (4.1), pj^ = P^ ± iP^\ and of the 
odd operators, 

Si = ^(F 1 +iF 3 ), S 2 = \ (F 2 + iF 4 ) , Q 1 = ±(F 7 + iF 5 ), Q 2 = ^(F 8 + iF 6 ), (4.6) 



P (±) p (±) 
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Table 7: Fermion-fermion anticommutation relations. Here S 1(2 ) = 2 U +> + (-)4 
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Table 8: Boson-fermion commutation relations 



and denote S\ t2 = S\ 2 , Qi j2 = Qi^ In terms of these complex combinations, the 
(anti) commutation relations including fermionic integrals are presented in Tables 7 and 8. 

From the (anti)commutation relations it follows that the operators and P±~^ act invariantly 
(in the sense of commutator) on every of the spin one-half multiplets (Qi,«?i), (Si, Qi), (Q 2 ,S 2 ), 
(S 2 , Q 2 ), where the first and second components of every multiplet are the eigenstates of 

with eigenvalues +^ and —\. Analogous multiplets for the bosonic operators J3 \ and P± ^ are 
(Qi,<5i), (<Si,Qi), (Q 2 ,S 2 ) and (S 2 ,Q 2 ). The eigenstates of the u(l) generator are the linear 
combinations of (4.6) and of the Hermitian conjugate operators: the operators (Si —S 2 ), (Si +S 2 ), 
(Qi — Q2) and (Qi + Q 2 ) are the eigenstates of the eigenvalue +1, while (Si + S 2 ), (Si — S 2 ), 
(Qi + Q2) and (Qi — Q 2 ) are the eigenstates of the eigenvalue —1. The anticommutators of the 
odd integrals produce linear combinations of H and the generators of the two chiral su(2) algebras 
modulo H. 



5 Supersymmetry reduced on the ground state 



As we have seen, the three possible gradings result in only two forms of the supersymmetry algebra 
given by the (anti)commutation relations with A = for T = R, and with A = 1 for V = 0-3 and 
T = Res. The reduction of these superalgebras on the subspace of scattering states with fixed 
energy E > fl 2 > effectively results, for both values of A, in the same superextension of the 
u(l) © su(2) © su(2) Lie algebra. However, the essential difference between the two forms reveals 
itself on the ground state of zero energy. 
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Indeed, with H = 0, the bosonic part of the superalgebra is reduced to the algebra u(l) © e(2) © 
e(2), where the first term corresponds to the integral T, and other two correspond to the two copies 
of the 2D Euclidean algebras generated by the rotation operator J% + \ the translation generators 

Pa + \ [Pi + \ P^] = 0, and by its analogs for the T = — 1 subspace. 

For A = 1, in accordance with Table 7, all the fermionic operators anticommute, and in this 
case we have also the commutation relations of the form [P, Q] = 0, [P,S] ~ Q, where we do not 
display the indices and do not distinguish fermionic operators from the Hermitian conjugate. 

On the other hand, for A = 0, the fermionic integrals commute with the translation generators, 
[P, S] = [P, Q] = 0, and anticommute between themselves in the following way: {5,5} = {Q, Q} = 
0, {S, Q} ~ P. Therefore, the difference between the grading T = Ron one hand, and the gradings 
r = 0-3 and T = R(T3 on the other hand, reveals itself on the ground state. 

It is interesting to note that the supersymmetry of the form similar to the present supersym- 
metry reduced on the ground state was analyzed by Gershun and Tkach [7] in the context of 
spontaneous supersymmetry breaking in 3+1 dimensions. The superalgebra used by them contains 
two parameters, a and b, (see Eq. (1) from [7]). The cases (a = 0, b 7^ 0) and (a 7^ 0, b = 0) 
correspond, respectively, to our cases A = and A = 1. 



6 Concluding remarks and outlook 



Although the system we have studied is simple, the existence of nontrivial (nonlocal) integrals of 
motion reveals a new super symmetric structure. A part of the nontrivial integrals of motion of the 
N = 2 superextended system comes from the hidden supersymmetry of the corresponding bosonic 
(spinless) quantum mechanical Dirac delta potential problem, while other conserved quantities ap- 
pear naturally with the addition of fermionic degrees of freedom. Combining both supersymmetries, 
the usual one generated by (2.3), and the hidden supersymmetry generated by (2.8), we find new, 
a priori unexpected nontrivial (nonlocal) integrals of motion (2.11). As a result, we enrich the 
knowledge on the system with a new property: its supersymmetric structure admits three different 
grading operators. With respect to these three different Z2-gradings, every of the three sets of 
nontrivial integrals (2.3), (2.8) and (2.11) (and other associated integrals obtained from them by 
multiplication with grading operators) is identified once as a set of bosonic generators, while two 
other sets are identified in this case as fermionic operators. 

Hidden supersymmetry exists also in a bosonic (spinless) Poschl- Teller (PT) quantum problem 
for integer values of the parameter of the system (for the details see ref. [5]), the appropriate limit 
of which corresponds to the Dirac potential problem. However, in PT problem, hidden N = 2 
supersymmetry is nonlinear [9, 10], of the order defined by the value of the parameter. Therefore, 
one should expect that such a hidden nonlinear supersymmetry has also to reveal itself in the 
simplest N = 2 superextended version of the PT system. However, due to the nonlinear character 
of the hidden supersymmetry, the resulting complete supersymmetry has to have a more rich nature 
in comparison with the Dirac delta potential case. Moreover, since the PT problem at those special 
values of the parameter (being reflectionless) carries some characteristics of a free particle, and in 
a certain limit produces a particle model of conformal mechanics [11], it would be interesting to 
investigate the PT system from the viewpoint of possible hidden superconformal symmetry. 

Another system that reveals hidden nonlinear supersymmery is the quantum periodic Lame 
problem [12], which reproduces the hidden supersymmetry of the PT system in a certain limit. 
Following the line presented in this work, we also are going to investigate the question on hidden 
supersymmetry in the superextended Lame system, which has some remarkable properties in the 
context of supersymmetry [13]. 

Let us note that the observed exotic supersymmetric structure with three different Z2-gradings 
can be effective in the search of hidden symmetries in some other quantum systems. In particular, 
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a general knowledge of this structure is helpful for identification of the unknown hidden symmetries 
in the broad class of the associated Lame systems [14], and in solving them. The key point there 
is that to every Lame system two different algebraisation schemes are related [15], with which two 
different sets of the non-diagonal integrals of the type (2.3) and (2.11) can be naturally associated. 
The detailed study of Lame and Poschl- Teller systems in the light of the revealed here nontrivial 
supersymmetric structure will be presented elsewhere [16]. 

Since the reduction of the observed hidden nonlinear sit(2|2) supersymmetry of the N = 2 
superxtended Dirac delta potential problem to the ground state of zero energy produces the 2D 
Euclidean analogs of two particular cases of a special superextension of the Poincare symmetry [7] , 
it would be interesting to investigate the origin of this latter symmetry. One could expect that it 
should appear in a special (contraction) limit of some superextension of the de Sitter symmetry. 

Finally, it would be interesting to investigate the question on existence of the hidden supersym- 
metry in the systems of many particles, in particular, with Dirac delta interactions, as well as in 
nonlinear integrable systems. 
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